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ob 1 Abstract 
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Motivated by the recent success of optical trapping of alkali bosons, we have 

O 

studied the zeroth sound modes of dilute Fermi gases with arbitrary spin-/, 
which are spin-S 1 excitations (0 < S < 2/). The dispersion of the mode (S) 
J> . depends on a single Landau parameter , which is related to the scattering 

lengths of the system through a simple formula. Measurement of (even a 
subset of) these modes in finite magnetic fields will enable one to determine 

00 . 

0^ . all the interaction parameters of the system. 
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Since the discovery of Bose-Einstein condensation (BEC) in dilute gases of alkali atoms 
], there have been experimental efforts to cool alkali fermions such 6 Li and 40 K down to 



the degenerate limit. In current experiments, alkali atoms are confined in magnetic traps, 



which confine only the spin states aligned with the local magnetic field. As a result, the 
spinor nature of the atom is suppressed. The recent success of optical trapping || , however, 
changes the situation. In optical traps, different spin components are degenerate in the 
absence of magnetic fields. One therefore has the opportunity to study dilute Bose gases 
with integer hyperfine spins (or simply "spin") / > and Fermi gases with spins / > 1/2. 
In a recent paper ||], we have discussed the structure of Cooper pairs of alkali fermions in 
optical traps. Since most alkali fermions have spin / > 1/2, their Cooper pairs can have 
even spin J ranging from to 2f — 1. The internal structures of these large spin Cooper 
pairs will give rise of to a great variety of superfluid phenomena. 

The purpose of this paper is to study a key normal state property of dilute fermi gases 
with general spin-/ in the degenerate limit, their collisionless or "zeroth" sound. We 
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shall show that in addition to the ordinary density mode, the system has additional modes 
corresponding to coherent inter-conversions of different spin species. These modes are the 
generalizations of the spin waves of spin 1/2 Fermi liquids. As we shall see, the dispersions 
of the zeroth sound modes contain the information of all the interaction parameters of the 
system, i.e. the set of s-wave scattering lengths {cLj} of two spin-/ atoms in the total spin J 
channel. Thus, observation of these modes will not only provide evidence of the degenerate 
nature of the system, but also information about the scattering lengths a s , and hence the 
existence of superffuid ground state as well as their transition temperatures. 

As in our previous study |4j], we shall focus on the homogenous case, i.e. without ex- 
ternal potential. This is a necessary step before studying trapped fermions. Moreover, it is 
conceivable that optical traps of the form of cylindrical boxes (rather than harmonic wells) 
be constructed in the future. In that case, the discussions here will be directly applicable. 
As in our previous work ||], our symmetry classification of the spin structure (which is a 
crucial step in our solution) also applies to arbitrary potentials. 

In addition to homogeneity, we shall also consider the weak magnetic field limit, i.e. 
when the Zeeman energy is much smaller than the kinetic energy of the system. These 
are the regimes where the spinor nature of the fermi gas is manifested most clearly. As 
demonstrated by the recent experiments at MIT H, this limit can be easily achieved by 
specifying the total spin of the system. Since the low energy dynamics of the system is 
spin conserving ||, the specified spin can not relax. The system therefore sees an effective 
magnetic field with which its spin would be in equilibrium, a field which can be much smaller 
than the external field B ext . In the following, we shall refer to this effective field simply as 
"magnetic field" B, with the understanding that it is a Lagrange multiplier that determines 
the total spin of the system ||. 

(A) Zero magnetic field : We begin with the linearized kinetic equation for the distri- 
bution function matrix 5n p in the collisionless regime 




(1) 
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Our notations in eq.(|I]) are the same as ref. 0. Here, n° is the T = Fermi function, 
v p = V p e p , 5h p is a (2/ + 1) x (2/ + 1) matrix in spin-space, [Sh p (r,t)] a p— J dxe~* p ' x < 
ipt{r — x/2, t)^ a (r + x/2, t) > where ip a is the field operator. The energy matrix Si describes 
the change in the Hamiltonian due to 5n, 

[feplap = J dr7a 7 ,/35(P,P')[^p']<57 ( 2 ) 

where dr' means dp'/(27r) 3 , and fa-y,/3s(p,p f ) ar e the Landau parameters which can be 
extracted from the Hamiltonian of the system derived by one of us ||. It is shown in 
ref. that only the lowest hyperfine states (with spin /) will remain in the optical trap 
and that the interactions between these spin-/ atoms are spin conserving, of the form 

Hint = I E 9j E [ drOj M (r)Oj M (r), (3) 

Z J=0,2,.. M=-J J 

where JM (r) = J2 a p < JM\ffaf3 > ipa^ippir), and < JM\f 1 f 2 a(3 > are the Clebsch- 
Gordan coefficients for forming a spin- J object from a spin-/! and a spin-/ 2 particle [Q, 
g r = A-Kf^aj/Mp, and Mp is the mass of the atom. Pauli principle implies that only even 
J's exist in eq. ([$]). Evaluating < H int > in Hartree-Fock approximation, and using the fact 
that 5 < if)+(r)i()p(r) >= J dT[5n p (r)]p a , it is straightforward to show that 

2/-1 J 

fay,3S = 2 E 9j E <JM\ff 1 axJM\ff5(3>, (4) 

J=0,2,.. M=-J 

which is momentum independent as a result of s-wave interaction. Note that if g 3 < 0, the 
system will have a superfluid instability towards spin- J Cooper pairs at a sufficiently low 
temperature Tj J ) 0. Our discussions for negative g/s therefore applies to temperatures 
above Tj J ) but low enough so that the Fermi gas is degenerate. Before proceeding, we 

— -gf-) v p , which turns eq.QT|) and (0) into 

d t vp + v p ■ V (z/ p + 5ep) = 0, 5e p = N F f ait p S < v v >, (5) 

where iV^ = mkp /2% 2 h 2 is the density of state of a single spin component at the Fermi 
surface, kp is the Fermi wavevector, and < (..) >= $-£-{■■) denotes the angular average 
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over the Fermi surface. Note that the quasi particle energy Se p is isotropic in k-space as a 
consequence of the s-wave interactions between the particles. 

Next, we note that a rotation 6 in spin space will cause a change a a — > D^}(9)ap : where 
D^p is the rotation matrix in the spin-/ space. This implies z> — > v' = 

£(/)££>(/)+. /From 

eq.(|2]) and (f|), one can see that 5e p transforms the same way, 5e p — > D^5e p D^ + . Since 
v is made up of two spin-/ objects, it can be decomposed into a sum of spin-S" quantities 
Z> (5,M) ) which transform as [i)U)p{8M)f){/)+] ap = [u( S ' M \ p D { j^ M , where < S < 2/, 

— S < M < S. The solution of this equation is easily seen to be [Sh^ s ' M ^] a p oc< fa\SfM(3 >. 
We then have the representation 



1 M*>t))aB ^ 



£ < fa\SfM/3 > 

S.M 



V [M r >*)] a/3 j 

Substituting eq.(0) into eq.(E|) and using the identity 



V 4 S ' M, (r,<)/ 



(6) 



£ < JM'\Sha X JM'lffS/3 >< fi\SfM~i > 

■ySM' 

= {-ff-\2J + l)W(ffff; JS) < fa\SfM(3 > (7) 
where W is the Racah coefficient ||, eq.(H) becomes diagonal in the (S,M) modes, 

fly>i fl ^+v p .vr^ + 4^) = o, (8) 



6e fM) = F (S) < V [S,M) >j (9) 

F {S) = - £ ^ i (2J+l)W(////;JS), (10) 

J=0,2,... 7F 

where we have used the fact that N F g } = ^kpajj^ and (— 1) 2 ^ J = —1 in obtaining eq. (PH). 
Eq.(||) and (g) imply that 

(| + v P ■ v) 4 5 ' M) + F<*>Vp • V < M) >= 0. (11) 

which is precisely the equation for the ordinary zeroth sound mode with only I = spin- 
symmetric Landau parameter F/ =0 non-zero @ and is given by F^ s \ 



The dispersion relations of modes described by eg. (pTTj) is well known 0. They are 

1 = C qVF% . (12) 
J-i 2 w - qv F x 

The properties of the modes depend crucially on the sign of the parameters F^ s \ When 
> 0, one has a well defined propagating mode. When —1 < F^ < 0, the zeroth 
sound mode is Landau damped. When F^ s ' < —1, the system is unstable against spin- S 
distortions. Because of the dilute limit, k F a << 1 and hence \F^ S '\ < 1, stability against 
spin-S 1 distortions is guaranteed. 

It is instructive to consider some special cases : 

(i) The density modes i>( 5=0 ) for fermions with arbitrary spin-f : Using the fact that 
W(ffff; JO) = {-lff' J /{2f + 1), we have 

^ } % f9f 4 , n E (2J+l)k Faj (13) 

nV + 1) j= ,2,.. 

In particular, if there are no superfluid instabilities in all angular momentum J channel, 
then _F( 5=0 ) > and the density mode will not be Landau damped. 

(ii) Fermions with spin-1/2,3/2 and 5/2 : For / = 1/2, eq.(|T3"D reduces to the well-known 
results F( 5= °) = -F^ 5=1 ) = N F g = 2k F a Q /nh g. For / = 3/2, using the tabulated values 
of the 6j-symbols [[TTJ to calculate the Racah coefficients, we find F^ S=QS) = k F (ao + 5a2)/7r, 
F (s=i) = F {S=3) = _k F (a + a 2 )/vr, and F^ 5=2 ) = k F (a - 3a 2 )/vr. Thus the S = 1 and 3 
modes are always degenerate, and the degeneracy between the 5 = and S = 2 modes are 
lifted only by the interaction in the J = 2 channel. If there are no superfluid instabilities 
in any J channel, i.e. both ao, a 2 > 0, then S = 1 and S = 3 modes are always Landau 
damped. 

For large /'s there are no obvious systematics except for the S = result noted above. 
Modes for different 5"s are typically not degenerate, barring accidental values of g/s. In 
the case of / = 5/2 such as 22 Na and 86 Rb, we obtain F (0) = (2a /3 + 10a 2 /3 + 6a 4 )k F /n, 
F<V = -(2a /3 + 46a 2 /21-6a 4 /7)A; F /7r, F® = {2a /3 + a 2 /3-3a i )k F /Tr, ^ (3) = -(2a /3- 
29a 2 /21 + 19a 4 /7)A; F /7r, F^ = (2a /3 - 5a 2 /3 - a±)k F /7r, F^ = -(2a Q /3 + 25a 2 /21 + 
ai/7)k F /ii. 



(iii) Undamped zeroth sound modes in the dilute limit : Zeroth sound modes (S, M) with 
F^ s > > are propagating J7|. Because of the dilute condition k F a 3 « 1, we have « 1. 
In this limit, eq.([L^) can be integrated to give 

^)( g ) = g^(l + 2e- 2 e- 2 / F(S) ). (14) 

Since F^ << 1, the exponential term in eq.(fL4]) will have little contributions. The frequen- 
cies of zeroth sound for all S are essentially given by qvp. As a result, it will be difficult 
to obtain information of the interaction parameters from zeroth sound frequencies in zero 
field. On the other hand, we shall see that even a small magnetic field will cause significant 
changes in the zeroth sound dispersions, which lead to many observable features and enable 
one to determine all the interaction parameters. 

(B) Weak magnetic fields: When B ^ 0, the kinetic equation eq.([I]) will have an 
additional term X p = f[e p ,n p ] on the left hand side 0. At the same time, the equi- 
librium distribution function and quasiparticle energy (denoted as h° B and e° B respec- 
tively) are altered from the zero field values (n° and e°). The difference 5e p = e pB — e° is 

IKU= -vBF*f3+JdT'fajA Sn U^ with 5 K= K,b -K= (IJ) K- These two relations 
imply 

[6e° p U = ~^BF z ap - N F f aJ AKh, (15) 

where /i is the magnetic moment of the atom. The solution of eq.(|15|) is 5e p = cF z . Using 
the fact that (F z ) a p oc< fa\lf0{3 >, it is easily to show from eq.(|6]), fllOD , and flT5|) that 

5e° p = -fiB cS F z , B cS = B/(l + F (1) ). (16) 

Linearizing about the equilibrium configuration h pB and e pB , we have X p = 
| (\Si p , 5hp] + [5e p , 5«- p ])- ^From the definition 5h p = — (-^p) %> the relation 5h° p = 
5e° p , and the property (a - f3)[v p ] a p = Esm < fa\SfM(3 > Mvf M) which follows 
from eq,(^J), we have 

KJoU -U^)T,< S-A8SMP > (SiM) + (17) 

\ Ut P/ SM 
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where fl = fiB cS /%. With this additional term on the left hand side of eq.(|]) and repeating 
the procedure as before, we have find that eq.(|9|) remains unchanged whereas eq.(H) becomes 

d t vf M) + [v p • V - iQM] + 4 5 ' M) ) = 0. (18) 

Thus, the zeroth sound modes can still be classified by the quantum numbers (S, M) in the 
weak field limit. The equation for the dispersion now becomes 

i = F (s) f 1 dx qvpx-nM 

J-i 2 uj + VLM - qv F x V ' 

which upon integration gives 

1 uj uj + VlM + qv F fiB 

- ln —FTT7 — - !» Q = i ; F77T ( 20 ) 



2qv F uj + flM -qv F ' 1 + F« 

Since the collective modes are excitations above the ground state, we only need to study 
the uj > solutions of eq.(^). In the following, we shall discuss only the zeroth sound 
modes that are not Landau damped, which requires \oj + QM\ > qv F in eq.(^). While 
many features of these propagating modes can be obtained analytically, we first display the 
numerical solutions of eq.([2t]) for S = 3/2 with > and F^ s ' < in fig.l and fig. 2 
respectively. The notable features of these modes are : 

(iv) Zeroth Sound modes near q = : Near q = 0, qv F /\QM\ « F^ s \ it is easily seen from 



eq.(g§ that |TT 



J s > M \ q ) = -QM{1 + F^) 



(21) 



Since uj( s < m *) > 0, only VtM < modes can be excited at q = 0. Note that all q = modes in 
finite field are not Landau damped irrespective of the sign of F^ s \ ^From eq. (pi]) , one can 
also see that all u/ 5 ' _ l M l) modes increase (decrease) as q 2 for F^ > and < 0. (See also 
fig.l and 2). 

(v) The F^ > case : For F^ > 0, zeroth sound modes with VLM > emerge from uj = 
when q > q M = QM/vf- (See fig.l). Expanding eq.(PQ) about (uj = 0, q = qu), one finds 
the dispersion in this neighborhood is 



uj - (qv F - SIM) = 2fiMe-( 1+1/F(S) ) (2QM/w) (22) 

Another simple feature one can derive from eq.(|20|) is that as q increases so that qvp » 
\SIM\, uj {SM \q) -> uj {s) (q) - SIM. The dispersions for all (S, M) modes become parallel to 
uj = qvF-, with all SIM < 0(> 0) modes shifted up (down) by an amount of |Mf2|. (See 
fig.l). It is also straightforward to show that zeroth sound modes with SIM 7^ lie above 
the particle-hole continuum of that particular M state, i.e. uj > —SIM + qvp. 
(vi) The < case : We find that the SIM < modes decrease monotonically as q 
increases, and vanish at q. M , = \SlM\/vp in a manner similar to eq. (|2"2"D . The entire mode M 
lies below the corresponding particle-hole continuum, i.e., uj < —SIM — qvF- Solving eq. (p20|) 
graphically, one can also see that there are no solutions with uj > when SIM > 0, implying 
the absence of zero sound modes with SIM > 0. 

Determination of scattering lengths : For scattering lengths |aj ~ lOOa^ where a# is 
the Bohr radius, a Fermi gas of density ~ 10 13 cm _1 will have kpOj ~ 10 _1 , which implies 
{F^} ~ 10 _1 . As mentioned in Part (A), it will be hard to determine the scattering 
lengths a j from the B = zero sound modes for these values of {F < - S ' ) } because their 
small contributions. On the other hand, in the presence of magnetic field, different zeroth 
sound modes (S, M) are separated. Since the interaction contributions to the zeroth sound 
frequency at q = and to the critical wavevector qM are of the form 1 + F^ instead of the 
essential singularity form in eq.(|nj), their contributions should be measurable for kpa f ~ 
10 _1 or even smaller. Note also that there are only (2/ + l)/2 scattering lengths [a j; (J = 
0, 2, .., 2f — 1)] whereas the number of zeroth sound modes in finite field is J2 2 J=o(2S + 1) = 
(2/ + l) 2 . Even some of these modes may not be excited, (as in the case of F^ < 0), there 
are still more the conditions on a T provided by the zeroth sound frequencies than the number 
of a j themselves. Thus, it is possible to determine the entire set of scattering lengths {a^} 
from the zero sound dispersions. 

This work is supported by NASA Grant NAG8-1441, NSF Grants DMR-9705295 and 
DMR-9807284. 
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Figure Captions : 

Figure 1 : The zeroth sound mode for = 0.5 > 0, / arbitrary, and f2 = +0.2. From 
upper to lower, the curves correspond to M — —2, —1, 0, 1, 2 respectively. Both u>, qvF, and 
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Q are plotted with arbitrary units. The vertical intercepts of curves M = — 2 and —1 are 
w ( 2 -2) = 2fi(l + F^) and cu^' 1 ^ = Q(l + F^ 2 )) respectively. The horizontal intercepts of 
the M = 1 and M = 2 curves are gif^p = and = 2f2 respectively. Zeroth sound 

modes for other S have different number of branches but behave similarly. 

Figure 2: The zeroth sound mode for the case = —0.2 < 0, / arbitrary. From upper 
to lower, the curves correspond to M — —2 and —1. The value of Q and the expressions 
for the vertical intercepts are identical to those in figure 1. The horizontal intercepts are 
q-iVF = £1 and q~iVp = 2fl respectively. 
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